#include <stddef.h>

#include "cctk.h"
#include "cctk_Arguments.h"
#include "cctk_Parameters.h"

void initial(CCTK_ARGUMENTS);

void initial (CCTK_ARGUMENTS) {
  int i, j, k;

  int istart = 0;
  int jstart = 0;
  int kstart = 0;

  iend = cctk_lsh[0]-1;
  jend = cctk_lsh[1]-1;
  kend = cctk_lsh[2]-1;

  for (k=kstart; k<kend; k++)
    for (j=jstart; j<jend; j++)
      for (i=istart; i<iend; i++) {
        int v = CCTK_GFINDEX3D(cctkGH,i,j,k);
        // something = rho[v];
      }
}
